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We study the production of the lightest neutralinos in the radiative process e + e~ — » x?X?7 i n l° w 
energy supersymmetric models for the International Linear Collider energies. This includes the 
QO ' minimal supersymmetric standard model as well as its extension with an additional chiral Higgs 

singlet superfield, the nonminimal supersymmetric standard model. We compare and contrast the 
dependence of the signal cross section on the parameters of the neutralino sector of the minimal 
' and nonminimal supersymmetric standard model. We also consider the background to this process 

» I , coming from the Standard Model process e + e~ — > vuj, as well as from the radiative production of 

C> |i the scalar partners of the neutrinos (sneutrinos) e + e~ — > w*7, which can be a background to the 

' radiative neutralino production when the sneutrinos decay invisibly. In low energy supersymmetric 

models radiative production of the lightest neutralinos may be the only channel to study supersym- 
metric partners of the Standard Model particles at the first stage of a linear collider, since heavier 
neutralinos, charginos and sleptons may be too heavy to be pair-produced at a e + e~ machine with 
= 500 GeV. 
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INTRODUCTION 



Supersymmetry (SUSY) is at present the only known framework [l[ in which the Higgs sector of the Standard 
Model (SM), so essential for its internal consistency, is technically natural 0]. Supersymmetry is, however, 
' not an exact symmetry in nature. The precise manner in which SUSY is broken is not known at present. 
However, the necessary SUSY breaking can be introduced through soft supersymmetry breaking terms that 
£NJ ' do not reintroduce quadratic divergences in the Higgs mass, and thereby do not disturb the stability of the 
hierarchy between the weak scale and the large grand unified (GUT) scale. Such terms can typically arise in 
supergravity theories, in which local supersymmetry is spontaneously broken in a hidden sector, and is then 
| transmitted to the visible sector via gravitational interactions. A particularly attractive implementation of the 
idea of supersymmetry, with soft supersymmetry breaking terms generated by supergravity, is the Minimal 
Supersymmetric Standard Model (MSSM) obtained by simply introducing the supersymmetric partners of the 
of the SM states, and introducing an additional Higgs doublet, with opposite hypercharge to that of the SM 
■ Higgs doublet, in order to cancel the gauge anomalies and generate masses for all the fermions of the Standard 
Model [H, El- In order for broken supersymmetry to be effective in protecting the weak scale against large 
radiative corrections, the supersymmetric partners of the SM particles should have masses of the order of a few 
hundred GeV. Their discovery is one of the main goals of present and future accelerators. In particular, a e + e~ 
linear collider with a high luminosity C = 500/6 -1 , and a center-of-mass energy of = 500 GeV in the first 
stage, will be an important tool in determining the parameters of the low energy supersymmetric model with 
a high precision [1, H, 0, fl Hi • Furthermore, polarisation of the electron (and positron) beam can enhance the 
capability of such a linear collider [ic| in unravelling the structure of the underlying supersymmetric model. 

In the minimal supersymmetric standard model the fcrmionic partners of the two Higgs doublets {H\,H2) 
mix with the fermionic partners of the gauge bosons to produce four neutralino states ^9, i = 1,2,3,4, and 
two chargino states \f, j — 1, 2. In the MSSM with i?-parity (R p ) conservation, the lightest neutralino state is 
expected to be the lightest supersymmetric particle (LSP). The neutralino states of the minimal supersymmetric 
standard model with R p conservation have been studied in great detail, because the lightest neutralino, being 
the LSP, is the end product of any process involving supersymmetric particle in the final state. 



* Electronic address: 
t Electronic address: 

* Electronic address: 

§ permanent address 



rahulOimsc .res . in 
ppandita@nehu . ac . in 
sharma@imsc . res . in 



2 



However, the MSSM suffers from the so-called jx problem associated with the bilinear term connecting the two 
Higgs doublet superfields H\ and H2 in the superpotential. An elegant solution to this problem is to postulate 
the existence of a chiral electroweak gauge singlet supcrficld S, and couple it to the two Higgs doublet superfields 
Hi and H2 via a dimensionless trilinear term XH1H2S in the superpotential. When the scalar component of 
the singlet superfield S obtains a vacuum expectation value, a bilinear term XH1H2 < S > involving the two 
Higgs doublets is naturally generated. Furthermore, when this scalar component of the chiral singlet superfield 
S acquires a vacuum expectation value of the order of the SU(2)l x U(1)y breaking scale, it gives rise to 
an effective value of fi (fi e ff = X < S >= Xx) of the order of the electroweak scale. However, the inclusion 
of the singlet superfield leads to additional trilinear superp otential coupling (k/3)S 3 in the model, the so 
called nonminimal, or next-to-minimal [ill, H2T [13, [lil. [Ill [la. [Tl| , supersymmetric standard model (NMSSM). 
The absence of H1H2 term, and the absence of tadpole and mass couplings, S and S 2 in the NMSSM is made 
natural by postulating a suitable discrete symmetry [l8|, [l9[ . The NMSSM is attractive on account of the simple 
resolution it offers to the /1 problem, and of the scale invariance of its classical action in the supersymmetric 
limit. Since no dimensional supersymmetric parameters are present in the superpotential of NMSSM, it is the 
simplest supersymmetric extension of the Standard Model in which the electroweak scale originates from the 
supersymmetry breaking scale only. Its enlarged Higgs sector may help in relaxing the fine-tuning and little 
hierarchy problems of the MSSM [2(|, thereby opening new perspectives for the Higgs boson searches at high 
energy colliders [2ll . |22| , and for dark matter searches [23| . In the nonminimal supersymmetric standard model 
the mixing of fermionic partners of Higgs and gauge bosons [24], HH, [2(| produces five neutralino states Xi> 
i = 1,2,3,4,5, and two chargino states xf 1 j — 1, 2. Furthermore, because of the presence of the fermionic 
partner of the singlet Higgs boson, the neutralino states can have an admixture of this SU(2)l x U(1)y singlet 
fermion, thereby affecting the phenomenology of the neutralinos in the nonminimal supersymmetric standard 
model. 

The lightest neutralino state (£?) of MSSM or NMSSM, being typically the LSP, is stable and therefore, a 
possible dark matter candidate [27l . |28| . Since the neutralinos are among the lightest particles in low energy 
supersymmetric models, they are expected to be the first states to be produced at the colliding beam experi- 
ments. At an electron-positron collider, such as the International Linear Collider (ILC), they can be directly 
produced in pairs 

e + + e-^x°+xl (1.1) 

which proceeds via Z boson and selectron exchange [5^, [3(| • In collider experiments the LSP escapes detection 
such that the direct production of the lightest neutralino pair 

e + + e-^x°i+xi (1-2) 
is invisible. Therefore, one must look for the signature of neutralinos in the radiative process 

e + +e- ^x°i+Xi+l (L3) 

where the final photon is radiated off of the incoming beams or the exchanged selectrons. We note that this 
process is suppressed by the square of the electromagnetic coupling. However, it might be the first proceess 
where the lightest supersymmetric states could be observed at colliders. The signal of the radiative process 
p.3[) is a single high energy photon with the missing energy carried away by the neutralinos. T he p rocess (|I.3[) 
has been studied in detail in the minimal supersymmetric model [U [H, [33], H3, [35], [H, H3, [H, 0, 53, 53 • In 
these studies, different approximations have been used in calculating the cross section for (|I.3|) . and the focus 
has been on LEP ener gies and special neutralino mixing, especially the case where the neutralino is a pure 

photino 0, m m m eh, [33, in . 

More recentl y, c alculations have been carried out in the context of MSSM assuming general neutralino 
mixing [39], lifl |4l|. Some of these studies underline the importance of longitudinal [3^] and even transverse 
beam polarisations. On the other hand, the signature "photon plus missing energy" has been studied in detail 
by different LEP collaborations ALEPH j42j, DELPHI [H[, L3 [3, and OPAL [H, El- We recall that in 
the SM, the radiative neutralino process e + e~ — > z/i/7 is the leading process with this signature, for which the 
cross section depends on the number N v of light neutrino species [47| . This signature has, thus, been used 
to measure N v , which has been found to be consistent with three. Furthermore, the LEP collaborations have 
found no deviations from the SM prediction, and, therefore, only bounds on the masses of supersymmetric 
particles have been set [42|, 0, 0, 46] . This process is also important in determining collider bounds on a very 
light neutralino [48| . For a review of the experimental situation, see Ref. [4!| . 

Most of the theoretical studies on radiative neutralino production in the literature have been carried out in 
the framework of the minimal supersymmetric standard model. This includes calculations relevant to ILC with 
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a high ccntcr-of-mass energy, high luminosity and logitudinally polarized beams, as well as study of the SM 
background from the radiative neutrino production 

e + e~^v + D + j, (1.4) 

and the MSSM background from radiative sneutrino production 

e + e~ — > £> + i>* + -7. (1-5) 

It has been pointed out that the discovery potential of the ILC might be significantly extended if both 
beams are polarized, especially if other SUSY states like heavier ncutralino, chargino or even slepton pairs are 
too heavy to be produced at the first stage of the ILC at sfs = 500 GeV. 

In this paper we shall consider the radiative process (|I.3|) in the nonminimal supersymmetric standard model. 
As discussed above, the nonminimal supersymmetric standard model is an attractive alternative to the MSSM, 
which solves the /1 problem of the minimal supersymmetric standard model, and in which the weak scale 
originates from the supersymmetry breaking scale only. Furthermore, the NMSSM has five neutralino states 
in its spectrum, and there is an admixture of a singlet state in the neutralino states, which may affect the 
radiative neutralino production process (|I.3j) . On the other hand, the background processes from the SM and 
supersymmetry are not affected by the spectrum of neutralino states. We shall compare and contrast the signal 
for the radiative neutralino process in the NMSSM with that in the MSSM, and study in detail the dependence 
of the cross sections on the parameters of the ncutralino sector. 

The plan of the paper is as follows. In Sec. [Til we discuss the cross section for the signal process (|I.3[) in 
NMSSM, and compare and contrast it with the corresponding cross section in the MSSM. Here we recall the 
basic features of the neutralino mixing matrix in NMSSM and MSSM, and the couplings of the neutralinos 
relevant for our calculations, as well as the cross section and the phase space for the signal process. Here we 
also describe the cuts on the photon angle and energy that are used to regularise the infrared and collincar 
divergences in the tree level cross section. We then describe the typical set of input parameters for the NMSSM 
that are used in our numerical evaluation of the cross sections. The set of parameters that we use are obtained 
by imposing various experimental and theoretical constraints on the parameter space of NMSSM. On the other 
hand, for MSSM we use the typical benchmark parameter set of the Snowmass Points and Slopes la ( SPS la) 
scenario [5l|, except when otherwise indicated. We analyse numerically the dependence of the cross section on 
the parameters of the neutralino sector, and on the selectron masses. In Sec. IIIII we discuss the backgrounds 
to the radiative neutralino production process (|I.3|) from the SM and supersymmetric processes. Here we also 
define a statistical significance for measuring an excess of photons from radiative neutralino production over 
the backgrounds. We summarise our results and conclusions in Sec. HVl 



II. RADIATIVE NEUTRALINO PRODUCTION 

A. Neutralino Mass Matrix, Lagrangian and Couplings 

In order to calculate the cross section for the radiative production of neutralinos 

e-( Pl ) + e+( P2 ) -> xS(fci) + XiCfo) + 7(9), (H.l) 

where the symbols in the brackets denote the four momenta of the respective particles, we need to compute 
the couplings of the neutralinos to electrons and the scalar partners of electrons, the selectrons. These can be 
obtained from the neutralino mixing matrix. To obtain the neutralino mixing matrix for the MSSM, we recall 
that the neutralino mass matrix obtains contributions from part of the MSSM supcrpotcntial 

I^mssm = nHiH 2 , {11.2) 

where H\ and Hi are the two Higgs doublet chiral superfields, and [i is the supersymmetric Higgs(ino) parame- 
ter. In addition to the contribution from the superpotential, the ncutralino mass matrix receives contributions 
from the interactions between gauge and matter multiplets, as well as contributions from the soft supersym- 
metry breaking masses for the gauginos. Including all these contributions, the neutralino mass matrix, in the 
bino, wino, higgsino basis (— iX', — iX 3 , i/jjj , ip% 2 ) can be written as [52I |53|| 

Mi —mz sin 9 W cos (3 mz sin 9 W sin (3\ 

M 2 mz cos 9 W cos (3 -m z cos 9 W sin j3 I m „.. 

Mussm — I _ m<7 s j n o_ ms ft ; -u [lL6) 



-mz sin 9 W cos (3 mz cos 9 W cos f3 
mz sin 9 W sin (3 —mz cos 9 W sin (3 



-I' 
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where Mi and M 2 are the U(1)y and the SU(2)l soft supcrsymmetry breaking gaugino mass parameters, 
respectively, and tan (3 = v 2 jv\ is the ratio of the vacuum expectation values of the neutral components of 
the two Higgs doublet fields Hi and H 2 , respectively. Furthermore, mz is the Z boson mass, and W is the 
weak mixing angle. In the CP conserving case, M is a real symmetric matrix and can be diagonalised by 
an orthogonal matrix. Since at least one eigenvalue of M is negative, we can use a unitary matrix N, the 
neutralino mixing matrix, to get a positive semidefinite diagonal matrix [52| with the neutralino masses m x q : 

N*M MSSM N- 1 = diag(m x o, to % o, m x o, m x o) . (11 A) 

We note that that the transformation Eq. (|II.4[) is only a similarity transformation if N is real. 
For the NMSSM, the relevant part of the supcrpotcntial is 

WWsm = \SHiH 2 -^S\ (II.5) 

where S is the Higgs singlet chiral superfield. In the basis ( — zA', — iA 3 , , i/ 1 ?/ >ips), the neutralino mass 
matrix for the NMSSM can then be written as [H, HI 



M, 



NMSSM = 



/ 


Mi 









— mz sin 6 W 


cos/3 


mz sin 6 W 


sin/3 


° \ 







M 2 




m z cos 6 W 


cos/3 


—mz cos 6>„, 


sin (3 







-mz sin 9 W 


cos /3 


mz cos 6 


wCos/3 







—Ax 




-Aw 2 


I 


mz sin 6 W 


sin/3 


—mz cos 6 


w sin/3 


—Ax 









— Aui 












-Av 2 




— Xvi 




2kx J 



(II.6) 



where < S >= x is the vacuum expectation value of the singlet Higgs field. As in the case of MSSM, we can 
use a unitary matrix N' to get a positive semidefinite diagonal matrix with the neutralino masses m x o 0, [25| : 

A^Mnmssm^'" 1 = diag(m x o, m x o, m x o, m x o m x g) . (11.7) 
The Lagrangian for the interaction of neutralinos, electrons, selectrons and Z bosons for MSSM is given by [52j 

C = (-^7T N n)fe p LxleR + ^ „ (^12 + t a n6 w N n )f e P RX ie L 
cos W V2 sin ^ 

+eZ Al / e7 ^ [-— — (\ - sin 2 9 U ) P L - tan 6 W P R ] f e + h.c, (II.8) 

sin fu, cos c/u, \Z ) 

with the electron, selectron, neutralino and Z boson fields denoted by / e , ei^, x®, and respectively, and 
Pr.l = h (l i 7' 5 ). The corresponding interaction Lagrangian for NMSSM is obtained from (|II.8p by replacing 
iVjj with Nlj . The different vertices following from (|II.8|) arc shown in Table |TJ 

The couplings of the lightest neutralino to electrons, selectrons and Z boson arc determined by the corre- 
sponding elements of the neutralino mixing matrix (iVy or N-j). For numerical calculation of the radiative 
neutralino cross section in the MSSM, we have chosen to work with the parameters in the SPS la scenario [5lj . 
The parameters of the SPS la scenario are summarised in Table HT1 However, since in the SPS la scenario the 
value of the parameters fi and M2 are fixed, we shall use a different set of parameters to study the dependence 
of the neutralino mass and the radiative neutralino production cross section on /1 and M% , and on the selectron 
masses. This set of parameters is shown in Table |IIH We shall call this set of parameters as the MSSM 
electroweak symmetry breaking scenario (EWSB) [541 ]. 

On the other hand, for the NMSSM we use a set of parameters that is obtained by imposing theoretical and 
experimental constraints on the parameter space of the NMSSM. The parameters that enter the neutralino mass 
matrix of the NMSSM are, apart from Mi and M 2 , tan /3, fj, (= A < S >= Ax), A and K. For Mi,M 2 and M 3 
we use the values which are consistent with the usual GUT relation M\/ai = M2I0L2 = M-^ja^- We note that 
for the MSSM in SPS la scenario, the value of the parameter tan j3 = 10. In order to remain close to the SPS la 
scenario of MSSM, we have chosen for our numerical calculations in NMSSM tan (3 — 10, whereas the rest of the 
parameters arc chosen in such a way that the lightest Higgs boson mass, the lightest neutralino mass and the 
lightest chargino mass satisfy the present experimental lower limits. We have also imposed on the parameter 
space of NMSSM the theoretical constraint that there is no charge and color breaking global minimum of the 
scalar potential, and that a Landau pole does not develop below the grand unified scale (Mqut ~ 10 16 GeV). 
For instance, for the parameter values tan/3 = 10, ^1 = 130 GeV, A = 0.55 and K = 0.44, we obtain the lightest 
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TABLE I: Vertices corresponding to various terms in the interaction Lagrangian (|II.8|) for MSSM. In addition we have 
also shown the vertices for selectron-photon and electron-photon interactions. The vertices for the NMSSM are obtained 
by replacing JVy with N^. 



Vertex Vertex Factor 

right selectron - electron - neutralino ~ te j JVfi Pl 

° cos W 11 

left selectron - electron - neutralino m le „ (N12 + tan 8 w Nh)Pr 

V2 sin W x ' 

neutralino - Z° - neutralino 4sin fl ^ e cos e ^ (\N l:i \ 2 - \N 14 \ 2 ) 7 M 7 5 

electron - Z° - electron ie-y^ [ sin9m 1 cos 9m {\ - sin 2 6 W ) P L - tan 9 w Pr] 

selectron - photon - selectron ie(pi + P2) 11 

electron - photon - electron ie^ 1 



Higgs boson mass of about 123 GeV, and the lightest neutralino boson mass of about 77 GcV, both of which 
are phenomenologically acceptable values. As another example, choosing tan/3 = 10 and fj, = 225 GeV results 
in the lightest Higgs boson mass of about 114 GeV, the lightest neutralino mass of about 132 GeV but with 
a charge and color breaking global minimum. Other choices give a Landau Pole below Mqut etc. for some 
particular set of free parameters. It is important to note that choice of parameters away from tan/3 = 10 is 
also possible, for example, tan/3 = 3, [i — 180 GcV, A = 0.3 and re = 0.1 produces a lightest Higgs boson mass 
of about 114 GcV, and the lightest neutralino mass of about 103 GeV. We have taken tan/3 = 10 for NMSSM 
simply because we wish to compare our results with MSSM for the typical SPS la scenario. The consequence 
of imposing these constraints on the parameter space of NMSSM, and the resulting masses for various particles 
for a particular choice of input parameters is summarized in Table IIVI 

Since the neutralino mass matrix depends on the parameters A and re, it is useful to study the possible 
values of these parameters, with all other parameters fixed, which satisfy the phenomenological and theoretical 
constraints discussed above. In Fig. [1] we show a plot of A versus re, with all other input parameters fixed as 
in Table IIVI and with the lightest neutralino, the lightest Higgs boson, and the lightest chargino mass as in 
Table ITVl with a variation of less than 5%. Fig. Q] shows the range of A and re values that are consistent with 
all the constraints discussed above for the set of input parameters in Table IIVI 

We note that for the set of input values in Table [TVl values of A < 0.4, with re < 0.22, lead to an unphysical 



TABLE II: Input parameters and resulting masses for various states in the MSSM SPS la scenario. 



tan/3 = 10 


H = 358 GeV 


M 2 = 192 GeV 


m = 100 GeV 


m x o = 97 GeV 


m x ± = 180 GeV 


mg R = 136 GeV 


m H = 400 GeV 


m v o = 181 GeV 


m ± = 383 GeV 

x, 


ms L = 196 GeV 


m h = 111 GeV 



TABLE III: Input parameters and resulting masses of various states in MSSM EWSB scenario. 



tan/3 = 10 


jU = 130 GeV 


Mi = 150 GeV 


Mi = 300 GeV 


M 3 = 1000 GeV 


At = 3000 GeV 


A b = 3000 GeV 


A T = 1000 GeV 


m x o = 97 GeV 


m ± = 119 GeV 
x, 


mg R = 143 GeV 


mp e = 194 GeV 


m v o = -141 GeV 


m ± = 330 GeV 

x? 


me L = 204 GeV 


m h = 120 GeV 
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TABLE IV: Input parameters and resuling masses of various states in NMSSM. 



tan/3 = 10 


H = 130 GeV 


Mi = 150 GeV 


M 2 = 300 GeV 


M 3 = 1000 GeV 


A t = 3000 GeV 


A 6 = 3000 GeV 


A T = 1000 GeV 


A = 0.55 


k = 0.44 


A A = 880 GeV 


A K = 10 GeV 


m x o = 77 GeV 


m ± = 121 GeV 
x, 


m gR = 149 GeV 


mp e = 194 GeV 


m v o = -158 GeV 


m ± = 334 GeV 

x, 


mg L = 209 GeV 


m h = 122 GeV 



global minimum. On the other hand, values of A > 0.57, with k > 0.44, lead to a Landau pole below the GUT 
scale. Thus, the allowed values of A and k, for the given set of input parameters, and for the fixed masses of 
lightest ncutralino, the lightest Higgs boson, and the lightest chargino, as in Table IIVI lie in a narrow range 
0-4 < A < 0.57 for 0.22 < n < 0.44. For definitcness, we have chosen to work with the values of A = 0.55 and 
K = 0.44 in this paper. These values correspond to the peak in the A versus k plot in Fig. [TJ 



X VS K 



0.45 



y 0.35 




0.25 



FIG. 1: Plot of A versus k for the set of input parameters in Table ITVl 



For the parameters of Table IIVI the composition of the lightest neutralino in NMSSM is given by 

N[j = (0.39, - 0.22, 0.57, - 0.59, 0.35). (II.9) 

From the composition (|II.9[) . we see that the lightest neutralino has a sizable singlet component, thereby 
changing the neutralino phenomenology in the NMSSM as compared to MSSM. For comparison, we also show 
the particle content of the lightest neutralino in MSSM 

N-Lj = (0.5, -0.21, 0.66, -0.51) (11.10) 

for the parameter set in Table HLT1 In Fig. [2] we have plotted the constant contour plots for the mass of lightest 
neutralino in NMSSM in the fi - Mi plane. We emphasize that the choice of \i and M2 values in this plot 
have been taken to be consistent with phenomenological and theoretical constraints as described above. For 
comparison, we have also plotted the corresponding contour plots for MSSM in Fig. [3] with parameters as in 
Table HU 
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B. Cross Section for the Signal Process 

In NMSSM, and in MSSM, the process proceeds at the tree level via t- and u-channel exchange of right 

and left selectrons e^L, and via Z boson exchange in the s channel. The corresponding Feynman diagrams 
are shown in Fig. [4j The differential cross section for can be written as [H, [55[ 

* = s^np^^o* +»-*■- (...ID 
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where p/ and Ef denote the final three-momenta ki, k2, q and the final energies E Xl , E X2 , and E 1 of the 
neutralinos and the photon, respectively. The squared matrix element \M.\ 2 in (jll.lip can be written as [33[ 

\M\ 2 = ]TT 4 „ (11.12) 

where TV, arc squared amplitudes corresponding to the Fcynman diagrams in Fig. [5] An average over intial 
spins and a sum over the spins of the outgoing neutralinos, as well as a sum over the polarizations of the 
outgoing photon is included in TV, . The phase space in is described in (33j . 



e (pi)- 



e+(P2)- 



wwvjM e (pi)- 

Xl(fcl) 

XI (fea) e+(p 2 )- 



— 4 K/vwv. 



-Xl(fel) e~(pi)- 
-XI (fe) 

»7(9) e+(p 2 )- 



-Xl(fcl 



en T 

-« : Xi (*a 



e (pi)- 



ww\-7(q) e (p x )- 

Xl(fc2) 



-Xl(fcl) e+(p 2 )- 



-Xl(fe) e (pi)- 



Xi(fci) 

*ww\. 7 (<j) e + (p 2 )" 



-Xl(fc 2 ) 



-« 1 Xl(fel) 




e (pi)- 



-Xl(*:i) e (pi)- 

-xi(fe 2 ) 



e + (P2) < k/vw\. 7 (q) e+(p 2 )- 



» : Xl(fel) e (pi)- 

Si t 

XiCte) e+(p 2 )- 



E JV\AAA. 7 (q) 

e" " 

Xi(fo) 

Sit 
1 Xl(fel) 



e (pi)- 



-Xl(k 2 ) e (pi)- 
-Xi(fel) 



e + (p 2 ) < 'vwwv-j(ij) e+(p 2 )- 



* 1 X1O2) 

] Xl(fcl) 



FIG. 4: Feynman diagrams contributing to the radiative neutralino production e + e 



X*txh- 



1. Numerical Results 



The tree- level cross section for radiative neutralino production (jll.lj) . and the background from radiative 
neutrino and sneutrino production, p.4[) and (|I.5|) . have been calculated using the program CalcHEP (54|. The 
tree level cross sections have infrared and collinear divergences, which need to be regularized. To do this we 
define the fraction of the beam energy carried by the photon as x = £? 7 /i?beam, where yfs = 2£'beam is the 
center of mass energy, and E 1 is the energy carried away by the photon. We then impose the following cuts 
on a;, and on the scattering angle (9 7 of the photon [33l. [Hoj: 

0.02 < x < 1 - p^ 3 -' ( IL13 ) 

beam 

-0.99 < cos 6> 7 < 0.99. (11.14) 
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FIG. 5: Photon energy distribution -777- for the radiative neutralino production for NMSSM (red solid line) and for 

dE- f 

MSSM (green dashed line) at s/s = 500 GeV. 



The lower cut on x in (|II.13j) corresponds to a photon energy E 7 = 5 GeV for -y/s = 500 GeV. The upper cut 
of (1 — w 2 o/£ ; boam) 011 x corresponds to the kinematical limit of radiative neutralino production. 

In order to implement the cuts on the photon energy in the calculation of the cross sections, we have taken 
the mass of the lightest neutralino in NMSSM to be m x o = 77 GeV for the parameter set shown in Table HVl 
For MSSM, we take m x o = 97 GeV from the SPSla scenario. Using Eq. (|II.13[) we get a fixed upper limit 
£™ = 226 GeV for NMSSM and £™ ax = 212 GeV for MSSM at = 500 GeV for the photon energy We 
have used exclusively these cuts for both signal and background processes. Two different mechanisms have 
been chosen for the production of neutrinos in the background process (|I.4|) , one "with upper cut" and another 
"without upper cut," for obvious reasons. We note that at = 500 GeV and for m x a > 70 GeV, this cut 
reduces a substantial amount of the on-shell Z boson contribution to radiative neutrino production process. 



2. Photon Energy (E 1 ) Distribution and Total Beam Energy (^/s) Dependence 

Using the procedure described above, we have calculated the energy distribution of the photons from radiative 
neutralino production in NMSSM and in MSSM SPS la scenario, respectively. This is shown in Fig. [51 where 
we compare the energy distribution of the photons in the two models. We note that the mass of the lightest 
neutralino is larger in MSSM SPSla scenario than in NMSSM for the set of parameters in Tables ITT1 and WV\ 
respectively, and consequently the upper cut for MSSM is lower than that for NMSSM. We also show the total 
beam energy y/s dependence of the cross sections in Fig. [S] for NMSSM and MSSM SPS la, respectively. Due 
to a smaller value of the mass of neutralino in NMSSM compared to MSSM SPS1 a, the total cross section 
in NMSSM is less than that in MSSM SPS la. From the kinematical endpoint £™ ax = E™*(m x a) of the 
energy distribution of the photon from radiative neutralino production, the neutralino mass can, in principle, 
be determined for yfs = 500 GeV. We note that although the shape of the energy distribution and the total cross 
section arc similar in NMSS and MSSM, the numerical values for NMSSM are considerably smaller as compared 
to the MSSM. This is primarily due to the fact that the lightest neutralino in NMSSM has a significant singlet 
component, thereby reducing the cross section. 
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FIG. 6: Total energy y/s dependence of the cross sections a for radiative neutralino production e + e — > x?X?7 f° r 
NMSSM (red solid line) and for MSSM SPSla scenario (green dashed line). 
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FIG. 7: Total cross-section a for the radiative neutralino production versus fi for NMSSM (red solid line) and for MSSM 
in the EWSB scenario (green dashed) at y/s = 500 GeV. For NMSSM fi = X < S >= Xx. 



3. Dependence on /i and M2 



Since the neutralino mass matrix, and hence the lightest neutralino mass, depends on [i and M2, it is 
important to study the dependence of the radiative neutralino cross section on these parameters. In the 
nonminimal supersymmctric standard model, /i (= A < S >= Air) and M% are independent parameters. We 
have, therefore, studied the cross section c(e + e~ — > XiXil) as a function of fi and M2 independently. In Fig. [7] 



11 



Neutralino Production 



0.01 



0.001 



a. 



1e-04 



1e-05 



I 


I I I I 




NMSSM 




-^MSSM EWSB 


I 


i i i i 



300 



400 



500 



600 

M 2 [GeV] 



700 



800 



900 



FIG. 8: Total cross section a for the radiative neutralino production versus M2 for NMSSM (red solid line) and for 
MSSM in the EWSB scenario (green dashed) at yfs = 500 GeV. 

we show the \x dependence of the total cross section for the radiative production of neutralinos for NMSSM as 
well as MSSM (EWSB). We recall that in the MSSM SPS la scenario these parameters are fixed. As is seen 
from Fig. [3 the total cross section increases with /i. The plot of total cross-section versus [i in Fig. [T] is plotted 
in the range fj, € [120, 200] GeV in NMSSM and in MSSM (EWSB). Note that the parameter values are chosen 
so as to avoid color and charge breaking minima, absence of Landau pole, and the phcnomcnological constraints 
on different particle masses. Furthermore, in Fig. [5J we show the M2 dependence of the total cross section for 
radiative neutralino production for NMSSM and MSSM. The total cross-section decreases with increasing value 
of M%. The graph of total cross-section versus M2 in Fig. [8jis plotted for the interval M2 £ [300, 900] GeV in 
NMSSM and in MSSM (EWSB) so as to satisfy the theoretical and phenomenoloigcal constraints described 
above. It is clear from Fig.[SJthat the radiative neutralino production cross section in MSSM decreases sharply 
as compared to that in NMSSM. 

4- Dependence on the Selectron Masses 

The cross section for radiative neutralino production er(e + e _ — > XiXi7) proceeds mainly via right and left 
selectron cr^l exchange in the t and w-channels. In the NMSSM and MSSM (EWSB), the selectron masses 
are independent parameters. In Fig. [5] and Fig. [TO] we show the dependence of total cross section of radiative 
neutralino production on the left and right selectron masses. The cross section is not very sensitive to the 
selectron masses for both models. Furthermore, the total neutralino production cross section is smaller in 
NMSSM as compared to MSSM (EWSB) as a function of left as well as right selectron masses. 

5. Photon Energy (E^) Distribution for the Radiative Production of the Second Lightest Neutralino 

As discussed above, the cross section for the production of the lightest neutralino in NMSSM is rather small 
compared with the corresponding cross section for the lightest neutralino in the MSSM. It may, therefore, be 
useful to consider the radiative production of the second lightest neutralino in the NMSSM. For the parameter 
set of Table IIVI the composition of the second lightest neutralino in NMSSM is given by 



N' 23 



(-0.085, 0.1, 0.69, 0.68, 0.19). 



(11.15) 
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FIG. 9: Total cross section a for the radiative neutralino production versus mg L for NMSSM (red solid line) and for 
MSSM in the EWSB scenario (green dashed) at yfs = 500 GeV. 
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FIG. 10: Total cross section a for the radiative neutralino production versus mi s for NMSSM (red solid line) and for 
MSSM in the EWSB scenario (green dashed) at <Js = 500 GeV. 



We have calculated the photon energy distribution for the radiative production of the second lightest neutralino 
in NMSSM for the set of parameters shown in Table llVl This is shown in Fig. [IT] For comparison we have also 
shown the photon energy distribution for the radiative production of the lightest neutralino in MSSM for the 
SPS la scenario. We see that the cross section for the production of the second lightest neutralino in NMSSM 
is much smaller than the cross section for the lightest neutralino in MSSM. 
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FIG. 11: Photon energy distribution for the radiative production of second lightest neutralino f° r NMSSM (red 
solid line) and for the lightest neutralino for MSSM (green dashed line) at yfs = 500 GeV. 



III. BACKGROUND PROCESSES 



A. The Neutrino Background 



The major background to the radiative neutralino production comes from the SM radiative neutrino 

production process [13, S3, [H, H3, [EH 

e + + e~ — > vi + V£ + 7 , I = e, /x, r. (IH-1) 

In this process ^ e are produced via i-channel VF boson exchange, and v e ^^ T via s-channcl Z boson exchange. 
We show the corresponding Feynman diagrams in Fig. 1121 The background photon energy distribution j^- 

and -y/s dependence of the cross section a for radiative neutrino production e + e~ — > is the same for both 
NMSSM and MSSM. As shown in Fig.[13]the photon energy distribution from the radiative neutrino production 
peaks at _E 7 = (s — m|)/(2y / s) w 242 GeV because of the radiative Z production(y / s > mz). This photon 
background from radiative neutrino production can be reduced by imposing an upper cut on the photon energy 
i max = E™ x /E hcam = 1 - m^ /El eam , see Eq. (llLT3l) . which is the kinematical endpoint E™ ax w 226 GeV of 
the energy distribution of the photon from radiative neutralino production 

m 2 x o = \(s- 2VIE7 aX ) • (III.2) 

In order to achieve this, one would have to separate the signal and background processes. This would be possible 
if the neutralino is heavy enough, such that the endpoint is removed from the Z° peak of the background 
distribution. 

In Fig. [14] we show the ^fs dependence of the total radiative neutrino cross section. Without the upper cut on 
the photon energy a; max , the background cross section from radiative neutrino production e + e~ — > is&j (green 
points in Fig. I14[) is much larger than the corresponding cross section with the cut (solid red line). However 
when we impose the cut, the signal cross section from radiative neutralino production is only about one order of 
magnitude smaller than the background in the case of MSSM, but is nearly three orders of magnitude smaller 
than the background in the case of NMSSM. 
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250 
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B. The Supersymmetric Background 



Apart from the SM background coming from (|III.1[) . the radiative neutralino production (jll.ip has a back- 
ground coming from the supersymmetric sneutrino production process [40L l59j 



+ v\ + 7 , 



e,fi,r. 



(III.3) 



The lowest order Feynman diagrams contributing to the process (|III.3j) are shown in Fig. 1151 This process 
receives i-channel contributions via virtual charginos for i/e^-production, as well as s-channel contributions 
from Z boson exchange for v e ^. T v* ^ T -production. In Fig. 1 161 we show the photon energy distribution 

for radiative sneutrino production e + e~ — ► i>i>*j at = 500 GeV. The total cross section for the radiative 
sneutrino production is shown in Fig. 1171 
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FIG. 14: The total energy y^s dependence of the radiative neutrino cross section cr(e + e~ — ► vv~f) with upper cut (red 
line) on the photon energy _E 7 , and of the radiative neutrino cross section without upper cut (green points) on the 
photon energy E-,, see Eq. (|II.13|I . 



Radiative sneutrino production (|III.3[) can be a major supersymmetric background to neutralino production 
pi.ip if sneutrinos decay mainly invisibly, e.g. via v — > Xa_ v - This leads to so called "virtual LSP" scenario [io| . 
However, if kincmatically allowed, other visible decay channels like v — > xfi^ reduce the background rate 
from radiative sneutrino production. For example in the SPS la scenario [HI, [6(| of the MSSM we have 
BR(£ e x?i/ e ) = 85%. 

Furthermore, neutralino production e + e~ — > X1X2 followed by subsequent radiative neutralino decay (r3l| 
X2 — * Xil 1S a ^ so a potential background. However, significant branching ratios BR(x° —* X?7) > 10% are 
only obtained for small values of tan/3 < 5 and/or Mi ~ Mi [4lT l62l. l63l]. Thus, we neglect this background, 
detailed discussions of which can be found in Refs. [62l. Ir33l. I&H . 
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FIG. 16: The photon energy distribution -j^- for the radiative sneutrino production e + e~ — > vv*~j at <Js = 500 GeV. 
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FIG. 17: Total energy y/s dependence of the radiative sneutrino production cross section a(e + e~ — > uu*j). 

C. Theoretical Significance 
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We now consider the question as to whether an excess of photons from radiative neutralino production can 
be measured over the SM background photons coming from radiative neutrino production. To quantify the 
excess of photons from the signal over the SM background photons for a given integrated luminosity C, we 
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FIG. 18: Theoretical significance S = -j=^=\fll for the radiative neutralino production versus /i for NMSSM (red 
solid line) and for MSSM in the EWSB scenario (green dashed). The value of y/s = 500 GeV. 
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FIG. 19: Theoretical significance S = ^ a \ a ^, f° r the radiative neutralino production versus M2 for NMSSM (red 
solid line) and for MSSM in the EWSB scenario (green dashed). The value of = 500 GeV. 



consider the theoretical significance [5C 



S 



iVs 



y/N S + N B V^ + CTB 



(III.4) 



where Ns = aC is the number of signal photons, and Nb = &b£ is the number of photons from the SM 
background of radiative neutrino production process. 
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The processes e + e~ — > XiXi7 ancl e + e~ — > depend significantly on the beam energy only near threshold 
in most of the parameter space for ^J~s = 500 GeV and C = 5 x 10 5 pb _1 . In Fig. [TH1 we show the \i dependence 
of the theoretical significance S for NMSSM. For this plot we have chosen A\ = 0, and varied the parameter 
fx in the interval fx e [110,220] GeV. All other parameters for the NMSSM are chosen as in Table [TV] For 
comparison we have also also shown in Fig. [T5] the theoretical significance S for the MSSM with parameters 
as in Table |TTT] for the EWSB model. We note that in NMSSM for fj, < 110 GeV, and with other parameters 
as described above, the Landau pole develops below the GUT scale, and the lightest chargino and the lightest 
ncutralino masses arc below the experimental lower bounds. On the other hand, for fj, > 220 GeV, the NMSSM 
develops unphysical global minima. Wc note from this analysis that Significance of S ~ 4 can be attained for 
fi ~ 220 GeV. 

Wc have also studied the behavior of the theoretical significance as a function of the SU(2)l gaugino mass 
M 2 . In Fig. [HI we show the M 2 dependence of the theoretical significance S for NMSSM and MSSM in the 
interval M 2 £ [300, 900] GeV. As in the study of the \x dependence, for NMSSM, we have chosen A\ = 0, and all 
other parameters as in Table HVl Furthermore, for the NMSSM the interval for M 2 is chosen in order to satisfy 
the theoretical and experimental constraints. For M 2 < 300 GeV, the Landau pole for NMSSM develops below 
the GUT scale. On the other hand for values of M 2 > 900, there is a sfermion with negative mass squared 
in the spectrum of NMSSM. We note from Figs. [15] and that higher values of fi and lower values of M 2 
are favored for achieving higher values of the significance S. A theoretical significance of S = 1 would mean 
that the signal can be measured at a 68 % confidence level. On the other hand a significance of 5 is required 
for the detection of the signal. We note that a signal significance of about 4 can be achieved for NMSSM for 
values of /it ~ 220 GeV. Besides the theoretical significance, one must also consider the signal to background 
ratio Ns/Nb in order to judge the reliability of the analysis. It will be necessary to do a detailed Monte Carlo 
analyses to predict the significance. However, this is beyond the scope of the present work. 

IV. SUMMARY AND CONCLUSIONS 

The nonminimal supersymmetric standard model is an attractive low energy supersymmetric model which 
solves the \i problem of the minimal supersymmetric standard model. We have carried out a detailed analysis 
of the radiative neutralino production e + e~ — * x?Xi7 m the NMSSM for the International Linear Collider 
energies and compared it with the corresponding results in the MSSM. This process has a signature of a high 
energy photon and missing energy. Wc have obtained a typical set of parameter values for the NMSSM by 
imposing theoretical and experimental constraints on the parameter space of NMSSM. For the set of parameter 
values that we obtain in this manner, the lightest neutralino in NMSSM has a significant admixture of the 
fermionic component of the singlet chiral superfield S. Using this parameter set, we have studied in detail 
the radiative neutralino production cross section in NMSSM for the ILC energies with unpolarized e + and 
e~ beams. For comparison, we have used the SPS la benchmark scenario for the MSSM. We have also 
calculated the background to this process from the SM process e + e~ — > vv"f, as well as the background from 
the supersymmetric process e + e~ —> vv*^j. All these processes have a signature of a highly energetic photon 
with missing energy. The photon energy distribution dcr/dEy, and the total cross section as a function of the 
total energy have been calculated for the NMSSM and for MSSM SPS la scenario at y 7 ! = 500 GeV using 
CalcHEP. Because of the admixture of a singlet in the lightest neutralino, the cross section as a function of 
energy for the radiative neutralino production is much lower in NMSSM than in MSSM. We have also studied 
the dependence of the cross section for radiative ncutralino production on the SU(2)l gaugino mass parameter 
M 2 and the Higgs(ino) mass parameter /x, as well as its dependence on the selectron (en, &l) masses in NMSSM, 
and compared it with the corresponding results in MSSM. In order to quantify whether an excess of signal 
photons, Ns, can be measured over the background photons, Nb, from radiative neutrino production, we have 
analysed the theoretical statistical significance S = No, / V As + Nb , and studied its dependence on M 2 and /i, 
the parameters that enter the neutralino mass matrix. The theoretical significance increases with the parameter 
fi for both NMSSM as well as MSSM, whereas it decreases as a function of M 2 . The decrease is especially sharp 
in the case of minimal supersymmetric standard model. It may be interesting to study whether the signal for 
radiative neutralino production in NMSSM can be enhanced by using polarized beams (65| . 
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